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This paper contains suplementary information for l| . The supplementary information is devoted 
to three main issues. In section I we describe the fabrication process; in section II we present the 
derivation of the Hamiltonian of the system and provide a more detailed discussion about the 
properties of the microbridges; in section III the hysteretic response of the resonator and the effect 
of heating are discussed. 

I. FABRICATION PROCESSES 

The fabrication process starts with a high resistivity Si substrate coated with SiN layers of 
thickness 100 nm on both sides. A 150 nm thick Nb layer is deposited on the wafer using magnetron 
DC sputtering. Then, e-beam lithography and a subsequent liftoff process are employed to pattern 
an Al mask, which defines the SSR and the SQUID leads. The device is then etched using electron 



cyclotron resonance system with CFa plasma. T 



400 Focus Ion Beam (FIB) system at accelerating voltage of 30 kV and Ga ions 

current of 9.7 pA. The outer dimensions of the bridges are about 150 x 50 nm. However, the 
actual dimensions of the weak-links are smaller, since the bombarding Ga ions penetrate into the 
Nb layer, and consequently, suppress superconductivity over a depth estimated between 30 nm to 



ie nanobridges are fabricated using FEI Strata 



50 nm 



m. 



II. DETAILED DERIVATION OF THE EFFECTIVE HAMILTONIAN 



The effective Hamiltonian of the closed system comprising the SSR and the SQUID 



3,Q is 



found using the hit me method thai was previously employed in Refs. lic| . ] fere however, we 
relax the assumption that the self inductance of the SQUID loop is small, and also the assumption 
that both junctions have the same critical currents. On the other hand, we assume that the 
inductance of the SQUID, which is denoted as L$, is much smaller than the total inductance of 



2 



the stripline Lt^t- This assumption can be justified by considering the fact that the measured 
angular resonance frequencies tu n of the first 3 modes (n 6 {1, 2, 3}) for all values of <& x (see Figs. 
2 and 3 in are very close to the values expected from a uniform resonator having length l^, 
namely ucot, where wx = tt/It^/LtCt- Moreover, the normalized flux-induced shift Auj n /nuJT in 
the angular resonance frequency of the first 3 modes is quite small and never exceeds 1CP 3 . Both 
observations indicate that the ratio Ls/Lt^t can indeed be considered as a small parameter. 

The resultant Hamiltonian of the closed system is given by TL = Wssr+^s (I), where TCssr is the 
SSR Hamiltonian and where Tis (I) is the SQUID Hamiltonian, which depends on the current / at 
the center of the SSR, namely, the current flowing through the SQUID. In terms of annihilation (A± 
and ^3) and creation (A\ and A^) operators for the first and third modes of the SSR respectively, 
the Hamiltonian TtssR can be expressed as 

W SS R = hu; T (JVi + 3iV 3 ) + U in , (1) 

where iVi = Ajvli and N3 = ^3^3 are number operators, 

V in = hs/^b 1 ? (e~ i ^ t A 1 + (2) 

represents the external driving, 7fi is the coupling constant between the 1st mode and the feedline, 
b™ is the amplitude of the driving pump tone, which is injected into the feedline to excite the first 
mode, and uj v is its angular frequency. 



A. The kinetic inductance of the nanobridges 



The Hamiltonian for the SQUID depends on the properties of the nanobridges. Due to the 
Ga ions implanted in the outer layer of the Niobium during the FIB process and the consequent 
suppression of superconductivity in that layer [(], ?|] , the weak links are treated as variable thickness 



nanobridges. 
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0, 



lavior of such a nanobridge is strongly dependent on the ratio Z/£ 
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12, 



18l |. where / is the bridge length and ^ is the coherence length of the Cooper 



pairs. The coherence length £ depends also on the temperature of the bridge. In the dirty limit £ is 
given by £(T) = 0.852-^/ £ // (Tc/T — 1) _1 where £0 is the size of the cooper pair and If is the 
mean free pathQ Q- The current-phase relation (CPR) of the bridges is periodic with respect 
to the gauge invariant phase 5 across the bridge. When l/£(T) <C 1, the nanobridge behaves like a 



regular Josephson junction (J J) with a sinusoidal CPR[2l[. However, as the ratio l/£(T) becomes 
larger, the CPR deviates from the sinosoidal form and can also become multivalued [l5|j. In case 



3 



the CPR is not multivalued the bridge can be approximately considered as a JJ having an extra 
kinetic inductance Lk • The effect of the kinetic inductance can be taken into account by replacing 
the screening parameter of the loop Pl = 27rA/ c /<J>o by an effective one given by 0l + A/3, where 
A/3 = 2irL K I c /<5> . 



In order to estimate A/3 we use Eqs. (47)-(49) and the data in Fig. 5 of Ref. 22|. For l/£ = 1.7 
the bridges' contribution is A/3 ~ 1 . As we will discuss below, both /3l and A/3 depend on the 
injected power P[ n that is used to excite the resonator due to a heating effect. However, for all 
values of Pi n that were used in our experiment, we estimate that the ratio A/3//3/, never exceeds 
the value 0.5 and thus the effect of kinetic inductance can be considered as small. Furthermore, 
the CPR remains a single valued function in the entire range of parameters that is explored in our 
experiments. Consequently, the nanobridges can be treated as regular JJs to a good approximation. 



B. The SQUID Hamiltonian 

In the following derivation we treat the nanobridges as regular JJs. We consider the case where 
the critical currents of both nanobridges are I c \ = I c (1 + a) and I C 2 = J c (1 — a) respectively, where 
the dimensionless parameter a characterizes the asymmetry in the SQUID. The Hamiltonian for 
the SQUID, which is expressed in terms of the two gauge invariant phases 5\ and 82 across both 
junctions, and their canonical conjugates p\ and P2, is given by 

2™ 2 p {pl+pl) 



+ E u(5 1 ,5 2 -J) 



(3) 



where uj p i = ^/I c /Cj^q is the plasma frequency, Eq = &qI c /it is the Josephson energy, and the 



dimensionless potential u is given by 



23] 



(1 + a) cos 5i + (1 — a) cos 82 



2 ~r $0 



(<5 a + 82) I C {h + 82 



4,1, 
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(4) 



where ( = $o/2/ c iT^T- 



C. Adiabatic approximation 



24j , the plasma frequency u> p i 



Due to the extremely small capacitance Cj of both nanobridges 
of the SQUID is estimated to exceed ITHz. Thus, the effect of the SQUID on the SSR, which 
has a much slower dynamics, can be treated using the adiabatic approximation 



25 



261 ] . Formally, 



treating the current / as a parameter (rather than a degree of freedom), the Hamiltonian Tig 
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can be diagonalized Hs\k(I)} = (/) \k (J)), where k = 0,1,2,..., and (k(I)\l(I)) = 5^. To 
lowest order in the adiabatic expansion the effective Hamiltonian governing the dynamics of the 
slow degrees of freedom correspon ding to the fast part of the system occupying the state \k (I)} is 



given by = Wssr + £fc (-0 [23, |28|. Furthermore, in the limit where the thermal energy k^T 
is much smaller than the typical energy spacing between different levels of H\ (~ fuv p i) one can 
assume that the SQUID remains in its current dependent ground state |0 (/)}■ For most cases this 
assumption is valid for our experimental parameters. It is important, however, to note that when 
the externally applied magnetic flux is close to a half- integer value (in units of <&o)> namely, when 
<f> x ~ (n + 1/2) $0) where n is integer, this approximation may break down. Near these points 
the potential u may have two different neighboring wells having similar depth. Consequently, near 
these points, the energy gap between the ground state and the first excited state can become much 
smaller than hw v \. On the other hand, the ratio between the height of the barrier separating the two 
wells (~ Eq) and the energy spacing between intra-well states (~ Huj p i) is typically Eq/Hlo p i ~ 100 
for our samples. Since the coupling between states localized in different wells depends exponentially 
on this ratio, we conclude that to a good approximation the inter-well coupling can be neglected. 
Moreover, in the same limit where Eo/hiOpi 1, one can approximate the ground state energy eo 
by the value of Equ at the bottom of the well where the system is localized. 

The current / at the center of the SSR can readably be expressed in terms of the annihilation 
and creation operators A\, A\ A 3 and A3. This allows expanding the current dependent ground 
state energy eq (I) as a power series of these operators. In the rotating wave approximation 
oscillating terms in such an expansion are neglected since their effect on the dynamics for a time 
scale much longer than a typical oscillation period is negligibly small. Moreover, constant terms 
in the Hamiltonian are disregarded since they only give rise to a global phase factor. In the 
present experiment the 1st SSR mode is externally driven, and we focus on the resultant dephasing 
induced on the 3rd mode. To that end we include in the effective Hamiltonian of the closed system 
in addition to the linear terms corresponding to the 1st and 3rd modes, also the Kerr nonlinearity 
term of the 1st mode, which is externally driven, and also the term representing intermode coupling 
between the 1st and the 3rd modes 

H eS = hu) X N x + hiu 3 N 3 + V in + KK x Nl + HXx, 3 NxN 3 . (5) 

The angular resonance frequency shift of the 1st and the 3rd modes, which is given by 

^1 ~ = ^3 ~ 3u^t = d 2 (eq/Eq) 
LOT ~ 3W T d {I/I c f ' 
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can be attributed to the inductance of the SQUID, which is proportional to the second derivative 
of £o with respect to /. On the other hand, the Kerr nonlinearity, which is given by 

Ki = C 2 f^i d 4 (fo/go) . , 

wi 2E d{I/I c f ' 1 ) 

and the intermode coupling, which is given by Ai^ = 9K±, can both be attributed to the nonlinear 
inductance of the SQUID [29], which is proportional to the fourth derivative of Eq with respect to 
I. 



D. Evaluation of u)\, u>3, K\ and A1.3 in the limit 0l -C 1 

The evaluation of the parameters U3, K\ and Ai,3 generally requires a numerical calculation. 

However, an analytical approximation can be employed when /3l <C 1. In this limit the phase 

difference 82 — 6% is strongly confined near the value 27r < I > :r /<I>o, as can be seen from Eq. This 

fact can be exploited to further simplify the dynamics by applying another adiabatic approximation, 

in which the phase difference 82 — 5\ is treated as a 'fast' variable and the phase average (5+ = 

{5 1 + 62) /2 as a 'slow' one. To lowest order in the adiabatic expansion one finds that for low 

frequencies to <C uj p i, namely in the region where the impedance associated with the capacitance 

of the JJs is much larger in absolute value in comparison with the impedance associated with the 

□ 

inductance, the SQUID behaves as a single JJ having critical current given by [301 ] 

is = 2i c ^l-(l-a 2 )sin 2 (tt$,/$o) • (8) 

Note that this approximation may break down when § x ~ (n + 1/2) $0 unless the asymmetry 
parameter a is sufficiently large. The relatively large value of a in our device (a ~ 0.5) ensures the 
validity of the above approximation. Using this result, it is straightforward to obtain the following 
analytical approximations: 



d 2 (eo/E Q ) 
d{I/Ic? 

d(I/Ic) 4 



ttI s 



3tt 2 I is 



(9a) 
(9b) 



which can be used to evaluate all the terms in Eq. (jSJ). 
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III. HYSTERETIC RESPONSE AND HEATING OF THE NANOBRIDGES 



As we discuss in the resonator exhibits hysteretic response to magnetic flux when the 
input power is relatively low. Such a behavior occurs, as can be seen from Eq. ([4j) above, when 
the screening parameter (3l is sufficiently large to give rise to metastability in the dimensionless 



potential u. A fitting of the model to the experimental data shown in Fig. 3(a) of [1| yields a 
value of 0l = 7.4. Another example of hysteretic response is shown in Fig. Q] below that shows 
data taken with another sample, which was fabricated using the same process that is described in 
the first section. The larger critical current in that sample yields a larger value of the screening 
parameter (3^ = 20. 

As is mentioned in the as the input power is increased the response becomes non- hysteretic. 
The gradual transition between the hysteretic region to the non-hysteretic one is seen in Fig. [2] 
below, which shows the difference in the measured resonance frequency of the first mode obtained 
from increased flux sweep (/ii nc ) and decreased flux sweep (/idee) at different input powers. Dark 
blue in the color map corresponds to no difference, namely to monostable regions, whereas in the 
red regions, where a large difference is observed, the system is bistable. As can be clearly seen 
from the figure, the bistable regions shrink as the input power is increased. The experimental 
results suggest that the critical current of the nanobridges drops as the input power is increased, 
and consequently the response becomes non-hysteretic due to the resultant smaller value of the 
screening parameter We hypothesize that the drop in the critical current occurs due to heating 
of the nanobridges by the input power. 

To estimate the effect of heating, we assume the case where the substrate is isothermal and that 



Forward Flux Sweep 



Backward Flux Sweep 




6.106 6.107 6.108 6.109 6.11 6.106 6.107 6.108 6.109 6.11 |S11| [dB] 



Frequency [GHz] 



Frequency [GHz] 



FIG. 1: Measured |5n| at input power P m = — 95dBm for forward (a) and backward (b) magnetic flux 
sweep. In this sample (3l — 20, and the response is highly hysteretic. 
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the heat is mainly dissipated down into the substrate rather than along the film [3l(. Moreover, 
we assume that most of the externally injected power into the resonator is dissipated near the 
nanobridges, where, the current density obtains its largest value. By estimating the heat transfer 
coefficient per unit area between each nanobridge and the substrate beneath it (100 nm SiN on top 
of high-resistivity Si) to be k ~ lWcm" 2 K -1 



32, 



331 ] and the area of the nanobridge to be A 



(50 nm) one finds that the expected temperature rise for P in = —70 dBm is AT = P\ n /An ~ 4K. 

Since heating is produced by AC current flowing through the nanobridges, it is important 
to estimate also the thermal rate, which characterizes the inverse of the typical time scale of 
thermalization, and is given by 7t = Ak/C, where the heat capacity C of the nanobridge is given 
by C = C v Ad, C v is the heat capacity per unit volume, and d is the thickness of the superconducting 
film. Using the estimate C v ~ 10~ 3 J cm -3 K _1 33] one finds 7t — 0.1 GHz. Since the frequency of 



the AC heating current is 1-2 orders of magnitude higher, we conclude that to a good approximation 
the temperature of the nanobridges can be considered as stationary in the steady state. 



-4 




WW l f „„ c - f ,dJ[ GHZ ] 

FIG. 2: The difference between the measured resonance frequencies obtained in the increasing flux sweep 
(/line) and the decreasing flux sweep (/idee) of the first (detector) mode. The dark blue areas correspond to 
monostablc regions, namely, the same resonance frequency is measured for both the increased and decreased 
sweep. The red indicates the regions where the system is bistable. 
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